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1. Introduction 

The BFSS |jl|] and IKKT [|| theories have generated recent interest in matrix models 
as tools for studying nonperturbative aspects of M-theory and type IIB superstring theory. 
In the BFSS formalism, M-theory formulated in the infinite momentum frame has been 
conjectured to be equivalent to the N ^ oo limit of a U{N) supersymmetric quantum 
mechanics written as a matrix model with N x N hermitian matrices. In the IKKT 
theory, type IIB Green Schwarz Superstring action in the Schild gauge has been related to 
a matrix model which has been derived as an effective theory of large N reduced model 
of 10 - dimensional Super- Yang Mills theory. In addition, there has been a recent surge of 
interest in noncommutative gauge theories after it has been realised that they appear as 
limits of D-brane open string theories in the presence of nonvanishing NSNS B fields , 

i- 

U (N) gauge theories in commutative spacetime can be studied, in the t'Hooft large N 
limit, by equivalent U{N) matrix models ^ where spacetime dependent fields are mapped 
to spacetime independent U{N) matrices to obtain a reduced model from the gauge theory 
and equality is obtained at the level of correlation functions of the field theory and the 
matrix model. The property of factorization of correlation functions of the gauge theory in 
the N —>■ oo limit, then guarantees equivalence of the gauge theory and the matrix model, 
provided the matrix model is either quenched (QEK model) or twisted (TEK model) 0. 
In the usual 't Hooft limit of the TEK model, a U{N) matrix model is equivalent to a 
U{N) field theory (in commutative spacetime) in a periodic box. In some recent papers 
0, however, it has been shown that there exists a limit of the twisted reduced models 
when they describe quantum field theory on noncommutative spacetime, in particular, 
a U{mn) twisted reduced model is equivalent to a U(m)e gauge theory (where 6 is the 
noncommutativity parameter) on noncommutative spacetime in the limit oi n oo. This 
limit is, clearly, different from the t'Hooft large limit ( in particular, nonplanar diagrams 
survive in this limit while they do not in the t'Hooft limit ) and, in a sense, it converts 
part of the noncommutativity of the matrix model into noncommutativity of spacetime 
coordinates, which are generated from the spacetime independent matrices in the reduced 
model by expanding a general matrix around a particular classical solution which exists 
only in the limit n ^ oo. Some other aspects of the noncommutative field theory and 
matrix model correspondence have beeen considered in 0. 

In this paper we use the relationship between noncommutative field theory and twisted 
reduced models to study Kosterlitz Thouless transition in noncommutative XY model. In 
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the commutative XY model, realised as a two dimensional system of interacting spins, 
there exist vortex configurations in which the order parameter is an integer multiple of 
the polar angle with the corresponding integer being the winding number of the vortex. 
Such vortex configurations are, however, not energetically favourable at all temperatures 
because of a competition between the energy required to form the vortex and the entropy 
of the vortex. Kosterlitz Thouless (KT) transiton 0] refers to the transition from a phase 
of no vortex to a phase of vortices and the temperature at which this transition occurs 
is called the KT transition temperature. This transition is characterstic of two spatial 
dimensions, since, in this case, both the vortex energy and the entropy of a single vortex 
scale as the logarithm of the system size. In noncommutative space one can, formally, 
construct the action for an XY model that reduces, in the appropriate limit, to the action 
for commutative XY model. Going to the reduced model one can, then, show that there 
exist vortex configurations for which the energy and entropy in two spatial dimensions scale 
as the logarithm of the system size. One can, by a calculation of the free energy similar to 
the commutative case, find an expression for the KT transition temperature, which turns 
out to be independent of the noncommutativity parameter 6. A brief comment on the 
nature of divergence of the noncommutative vortex action is, probably, in order here. The 
action scales as the sum of two terms one linear and the other a logarithm of the system 
size. By a suitable regularization, one can extract the logarithm piece for the action which 
gives the desired behaviour for the vortex energy. 

We organise the paper as follows. In section 2 we develop our notation and review the 
relation between noncommutative gauge theories and twisted reduced models. In section 
3 we discuss the noncommutative XY model and map it to the equivalent matrix model. 
In section 4 we find the vortex solution to the XY model and discuss Kosterlitz Thouless 
transition in the continuum limit of the model with infinite lattice size. In section 5 we 
study the same for the case of a finite lattice. Section 6 ends with some concluding remarks. 

2. Description of noncommutative field theory in terms of reduced model 

Field theories in noncommutative spaces can be studied by twisted reduced models 
obtained by a mapping of the fields on the noncommutative space to matrix valued op- 
erators. The usual matrix multiplication involving these operators then defines the way 
in which fields are multiplied in the corresponding noncommutative space. Integration of 
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fields over noncommutative coordinates is then equivalent to tracing over the correspond- 
ing matrix operators. On a noncommutative space of dimension D, the local coordinates 
are replaced by hermitian operators obeying the commutation relations 

[X'^,X^]=zr^ (2.1) 

where 6^^ are dimensionful real valied c-numbers. Fields are defined as functions on 
this coordinate space x^. The reduced model prescription is equivalent to finding matrix 
representation of the coordinate space, that respects the basic commutation relation ( ^TTI) , 
and of fields on this space. However, it turns out, that ( |2.1j ) does not have finite dimensional 
representations. For the purpose of doing explicit calculations one puts the field theory of 
interest on a lattice with some lattice spacing. Noncommutativity then implies that the 
lattice is periodic [|10[, [0, [0 ■ 



As an example, let us consider a scalar field defined on the noncommuting space 
with lattice spacing a and linear size L. We will specialize to two dimensions and hence 
/X = 1,2. We impose periodic boundary condition on the lattice and construct, formally, 
the matrices Fi = e*"^^ , F2 = e*"^^ ( where S = — ^ ), which generate translations 
by one lattice spacing in the directions X^ and X"^ respectively. In a periodic lattice it 
is possible to obtain finite dimensional representation of the F^ s even though the X^ s 
themselves are infinite dimensional. We consider an dimensional representation of the 
F^ s which satisfies the following relation, 

F1F2 = e^^FsFi (2.2) 

F;: = 1, (^ = 1,2) (2.3) 

where a = ^ . For an explicit construction of these matrices we refer to . All operators 
in this space can be written as linear combinations of integral powers of these matrices 
(since the set of all translations generate the full quantum space) 

N-l 

^= </'n.n,r;^^r^^e-^"^--"^ 5^^(fc.,fc,)e-^'=-^'' (2.4) 

where 4>n^ny = 4>{kxiky) and the k s belong to the discrete ( because the lattice size is 
finite) dual lattice. The corresponding field in noncommutative space is then defined as a 
fourier tranform of the scalar valued coefficient (t){k) 

<j>{x,y)=J2^ikx^ky)e-''^^'^ (2.5) 
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This gives a new coordinate space which one interprets as the semiclassical hmit of X^. 

We label the lattice, with Nq (L = Nqq) lattice sites along each dimension, by the 
set of dimensionless numbers {x,y) . The dimensional position coordinates on the lattice 
are then given hj x = xa , y = ya. The finite lattice size implies that the dimensional 
momenta and ky are quantized in inverse units of the lattice size, 

27112^: 2TTny . 

2; — 'IT? 5 — 'IT? (2-6) 

where and Uy are integers. We choose = A'^o and, hence, fourier expand the field (p in 
terms of dimensionless coordinates {x, y) and dimensional coordinates {x, y) on the lattice 
as follows, 

cPix, y)= '^n.n,e-^«("^^+-^^) = <^(x, y)=J2 ^(^- ^.)e-^^'^"+'-^^ (2.7) 

where (j)n^ny = 4>{kx, ky). 

To go to the continuum, we take the limit of vanishing lattice spacing while keeping 
the dimensional noncommutativity parameter 9 fixed. One can determine how 9 scales 
with the lattice spacing a and the lattice size L = Nqg = Na as follows. We consider the 
product of two matrix model operators (pi and (^2 given by 

N-l 

(2.8) 

Recalling the map from the reduced model to the noncommutative field theory we obtain 
the following equality 

9^'^{ki^k2y - kiyk2x) = a{nixn2y - niyn2x) (2.9) 
from which we obtain the following scaling equation for 9 

9^^ ^9=^ = ^ (2.10) 
2n 2n ^ ^ 



where we have used the relation L = Na. ( ^.10| ) is also verified from ( |2.1D and ( |2.2| ). 



The continuum limit is taken such that a ^ 0. Then, in order to keep 9 fixed, one has to 
take N ^ 00 with La = Na^ fixed. Thus, in the present case ( in section 5 we will take a 
different continuum limit such that the dimensional lattice size is fixed ), the continuum 
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limit forces one to go to the limit of infinite lattice size . In the continuum limit, the map 
between the noncommutative field (f) and the matrix model operator cf) remains the same 
as in ( |2.4[ ) and ( |2.5D with the sum over the discrete set k being changed to integration 
over continuous values of k. The derivative of the field 4> is then mapped to the following 
commutator in the reduced model 

id^^^[P^A] (2.11) 
where = {9~^) ^^^X'^ = B^^X^ and vector indices are raised or lowered with the metric 

Finally, we mention that the following relation holds between the integrals of fields on 
the noncommutative space and traces of operators in the reduced model, 

Tr4> = NcPoo = ^J^<PiS^^ y^^^ej y) (2-12) 

xy 



3. Non commutative XY model 

In this section we will consider the noncommutative XY model and try to understand 
the nature of Kosterlitz - Thouless transiton in it. The continuum action for the model is 
given by 

S = [ d^x{d^U+d^U) (3.1) 
9nc J 

where 

U*{x)*U{x) = l (3.2) 

Here U { U* is the complex conjugate of C/ ) is a scalar field on the two dimensional 
noncommutative space (x^, x^) and product of fields implies the star product. To see that 
this model is indeed the noncommutative analogue of the XY model one first uses ( |3.2|) to 
prove that U{x) can be written as U{x) = (e*^'-^'')*, where 6{x) is a new field. One can 
easily check that the action in (|3.1|) then takes the form 

1 r ?2 

s = / (fxidj *d''e + -(e* dj * d^e - dj * df'e * e) + oie^) + • ■ ■) (3.3) 

In the commutative limit, when star products are replaced by ordinary multiplication of 
functions, the action takes the form 5" = / (Px{dnO)^ which, one easily recognizes, 

9nc 
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is the form for the dimensionless energy in the corresponding statistical mechanical XY 
model with ${x) denoting the order parameter. 

The corresponding reduced model is obtained by replacing derivatives by commuta- 
tors, integrals by traces and star products by ordinary matrix products according to the 
prescription given above and is, therefore, given by the following action 

S = -^Tr[P^,U+][P^,U] (3.4) 

where 

U+U = UU+ = 1 and g^Na^ = g%c (3-5) 

We define creation and annihilation operators (we will see later that B is negative in our 
convention) as follows 

^ (Pi +iP2) and a+ = ^i=(Pi -iPa) (3.6) 



where 

[P^,P2] = -iB^2 = -iB (3.7) 

so that a and a"*" satisfy the usual commutation relation of the annihilation-creation oper- 
ators of a one dimensional harmonic oscillator 

[a,a+] = l (3.8) 

In terms of these variables the action takes the following form 

S ='^Tri[a+,U+][a,U] + [a,U+][a+,U]) (3.9) 
9 

and the equation of motion for the field U is given by 

[a+, U+[a, U]] + [a, U+[a+, U]] = (3.10) 

4. Vortex Solution in noncommutative XY model 

There exists a solution in this model which behaves as a vortex solution in the com- 
mutative limit. The solution is given by 

Uo = a^^^= (4.1) 
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This is the vortex solution for the higgs field in the abelian higgs model with an appropriate 



gauge field but also happens to be the vortex solution in our model in the absence of 
gauge field. To see that this is a vortex solution in the commutative limit we recall the 
definition of the creation and annihilation operators a and a"*" , 



a^ = ^-^{X^ + ^X') (4.2) 

Now maps to a new coordinate space interpreted as the semiclassical limit of X^. 
In this new coordinate space the solution Uq looks like 

Uo{z.z) = -—^ (4.3) 

[ZZj 2 

where z = [x^ — ix^) and z is the complex conjugate of z. Writing z in polar coordinates 
as z = rexp(z^), we regain the following behaviour of the above solution, namely, 

Uo{r,4>) = e'^^ (4.4) 

which is a vortex in two dimensions with winding number n. In the commutative XY model 
in two dimensions it can be easily shown, by a straightforward calculation, that the action 
for the vortex behaves as InL where L is the linear dimension of the two dimensional space . 
We evaluate the noncommutative action in the following way. We first consider a complete 
set of states of the number operator a+a and denote the individual states by |m >, with 
m running from to oo. We then evaluate the action as a trace over this complete set of 



states . However, it turns out that a complete analytic behavior of this action is not 
possible to obtain because of the infinite sum over the oscillator states. We regularize by 
introducing an upper cutoff R on the oscillator number m and evaluate the action as a 
sum only over states ranging from 1 > to | > . This prescription for regulating the action 
can be interpreted as an infrared regulation in the noncommutative theory because of the 
following reason. The spread ( root mean square deviation ), in postion space (x^,x^), 
of the wavefunction corresponding to the state \m > can be calculated to be VmO. The 
spread must be bounded by the system size L = y2TxN6. Equating the maximum spread 
with the system size we obtain R = 2nN . Since the spread depends on the oscillator 
number, imposing an infrared cutoff ( or fixing a particular system size ) implies a system 
size dependent cutoff on the harmonic oscillator level number m and is expressed by the 



relation L = \/R9 . Without the cutoff, we calculate the action for this vortex solution as 
follows, 

S = --^Tr[a^U^aUo + a^UoaU^ - {aa^ + a+a)] 



2^ OO CXD OO 

— [J^ y/m{m -n)+ ^ ^/m{m + n) - ^(2m+l)] 

^ m=n m=0 m=0 

— [2 ^ ^m(m + n) - ^ (2m + 1)] 

^ m=0 m=0 



(4.5) 



Here the limits on the first sum in the second line comes because 

OO OO 

< m\a^UQaUo\m >= m(m — n) (4-6) 

m=0 m=n 

and we have shifted the summed over integer to go from the second to the third line and 
used the following relations, 

a\m >= \^\m — 1 > , a'^\m >= y/m+ l\m + 1 > (4-7) 

Uo\m >= \m — n > , C^o"!™ >= \m + n > (4-8) 

Since we are interested in the way the action scales with the system size or, in other 

words its large distance behaviour, we shall try to investigate the behaviour, for large 
oscillator numbers, of the sum over oscillator levels in the action. 

Expanding the first term in the action as 

Vm(m + n) = m(l + i o^ + ---) (4-9) 

^ m 8 



we get 

9R °° 1 n2 

5 = f El("-i)-5^ + -l (^-lo) 

^ m=l 

The first term is linearly divergent and the second term is logarithmically divergent. 
One way to make sense of these sums is to use zeta function regularization, where 



OO ^ 

c(-) =j:^ (4-11) 



n—l 



While this is divergent for s = the analytic continuation of ({z) for arbitrary complex 
z is finite and equal to — |. Therefore, it is possible to assign a finite value to the first 



term. The second term is genuinely divergent and is in fact log divergent. Introducing the 
cutoff R on oscillator level m at this stage shows that the action diverges as In = In -g-. 
Thus we recover the usual behaviour that the action (or the dimensionless energy in the 
statistical mechanical viewpoint) for a vortex scales as the logarithm of the system size. 

One can follow a more rigorous method for finding the asymptotic expansion m. R = 
2nN of the action by using the Euler - MacLaurin formula, 

E/M = ^/(^)+/ dtfit) + C+[^f^'\R)-^f^'HR) + ^f('\R)----] (4.12) 

where C is a constant and B s are Bernoulli numbers. Using /(m) = 771(171 + n) in the 
above formula, we get for our action, 

2B rr^ 
S^-^ [R{n - 1) - —In R + 0{1/R) + ■ ■ ■] (4.13) 

This is the same expression as obtained from ( [4.10|) . 

A better regularization |jT^ is obtained by using the q analogue of the harmonic 
oscillator operators a and a+ denoted by A and respectively. 



where 

M _ -M 

[M], = i--^ (4.15) 

where q is 2i?-th root of unity and the usual oscillator limit is the one for g — > 1 or i? ^ 00. 
Replacing a, a"*" by A, respectively, and evaluating the sum over oscillators we get the 
following expression for the action, 

2B °° 

S = — E(2v/[m][m + n] - [m + 1] - [m]) (4.16) 

The infinite sum is automatically cutoff because [R]q = for q^^ = 1 and, hence, one 
does not have to put in a cutoff by hand. The action can, in principle, be evaluated, as a 
function of the cutoff R, by numerical methods. 

A small note on the Kosterlitz - Thouless (KT) transtion for this model is, probably, 
in order here. As noted earlier the field theory action of the vortex varies with R as 

S = -^InR = ^InR = ^In^ (4.17) 

ale ale 



The apparent negative sign in the first expression for the action may be surprising. How- 
ever, we note that with our choice of conventions 

r'^ = ^e'^^, Bf"" = {e-^y'' = -^^^'^ (4.18) 

from where it easily follows that S is a negative number 

5-5" = -^--^ (4-19) 

One can, of course, obtain a positive S by a change of convention which involves inter- 
changing the roles of X'^,X'^ but since the action ( |3.1| ) is symmetric in and x'^ it still 
remains positive. The entropy of a vortex depends on the number of available locations 
for the vortex. Because of the noncommutativity of the coordinates, a vortex cannot be 
localized in any region with area smaller than 9, since, the product of uncertainties of the 
two noncommuting coordinates is of order 9. This is true because from the relation 

[X\X^]=i9 (4.20) 

it follows that 

AX^AX^ > ^ (4.21) 

where AX^ denote the root mean square deviations of the coordinates X* taken over any 
state of the system. Therefore, for an area A = X^X"^ of order 9, the minimum uncertainty, 
A A, is of order 9 which implies that the centre of a vortex cannot be localized to a region 
of area less than order 9. Hence, the number of positions available to a vortex in a region 
of size is of order which implies that the entropy should scale as In^. Therefore, 
the dimensionless Helmholtz free energy is given, in terms of the couplings g and qnc by 

Hence, the couplings Qc and {gNc)c at which it is first favourable to produce vortices is 
given by 

9c = i9Nc)l = rrn^ (4.23) 

The critical temperature is obtained by using the following relationship between the field 
theory coupling gj^^j and the the statistical mechanical coupling J, namely. 



— (4-24) 
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where J is the strength of the spin-spin interaction that occurs in the expression for the 
spin-spin interaction energy E as follows 

E = J^cos{ei-ej) (4.25) 

where 9i refers to the order parameter at position i in a two dimensional lattice and the pair 
ij refers to nearest neighbours. With the correspondence ( [4.24|) the critical temperature 
is obtained as 

= JTrn^ (4.26) 



5. Effect of finite lattice size on KT transition 

In this section we will modify our previous calculations and take the continuum limit 
while keeping the lattice size fixed. This limit has been considered in fl^, fl^. We shall 
use the same dimensional representation of T ^ as in section 2. We expand the matrix 
model field ^ as follows, 

{N-l)/p 

<^= E </'n.n,rf''r?-e-i--^«^^^^^(/c.,fc,)e-^'='^^^ (5.1) 

where (pn^uy — (pikx^ky) and /3 = is an integer and A^o is the number of lattice sites 
along one dimension in the lattice. The upper limit on the sum above follows from the 
fact that = 1. For a lattice of spacing a and size L = Noa, the dimensional momenta 
kx and ky are again given as in ( p.6|) . The fourier expansion of the field (p in dimensionless 
and dimensional coordinates is then as follows 

(Pix.y) = J2 </'n.n,e-*"^('^^^+"''^) = (P{x,y) = Hk.ky)e-'(-^-^+''yy^ (5.2) 

f^x I'^y kx -jky 

The product of two matrix model operators (pi and (p2 is given by 

^1^2 = Yl ^n,.n,,<^n..n.,rr^rf^^e-^^^^^e-^('^--^^-— (5.3) 

where 

Nx = {nix + n2x)-, Ny = {niy + n2y) (5.4) 
11 



Recalling the map from the reduced model to the noncommutative field theory we obtain 
the following equation for 6 

The continuum limit ( a — > ) is taken such that the system size L = Nqg and the 
dimensional noncommutativity parameter 9 are kept constant. This implies the following 
scaling relations 

a=y^, (3=^V2^ (5.6) 

In this case, the following relation holds between the integrals of fields and traces of 
corresponding operators, 

N ^ ^ N f 

Trcf) = N(/)oo = -j^ 2^ (j){x, y) = — dxdy(f){x, y) (5.7) 

We now consider the noncommutative XY model action given by (|3.1| ). As usual, the 
corresponding reduced model is obtained by replacing derivatives by commutators, star 
products by matrix products and integrals by traces. Noticing that the map between 
derivatives and commutators remains the same as in ( p. 11 ), we get the following form for 
the reduced action 

5 = -lTr[P^,t/+][P^C/] (5.8) 

where = -^g^c- evaluation of the action along lines analogous to sections (3) and 
(4) then gives 

2B ^^'^ ^^'^ 
5 = — [2 ^ Vm(m + n) - ^ (2m + 1)] (5.9) 

^ m=0 m=0 

The upper limit on the sum follows from the fact that the spreads of the harmonic oscillator 
levels |m > ( which scale as \/m6 ) should be bounded by the size L of the system. 
Evaluating the sum in ( |5.9|) we obtain ( after subtracting the spurious linear dependence 
on LF' /6), for large LP' /6, the following expression for the action. 

The entropy is given by the same expression as in the infinite lattice case. The dimension- 
less Helmholtz free energy is, therefore, given by 

^- ^In— -In— = - In— -In— (5.11) 



2g^ e e 2Ngj,c 6 
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from which it follows that the critical values Qc and {gNc)c of the couplings above which 
it is favourable to produce vortices is given by 



9c 



{9nc)\ 



2N 



(5.12) 



Prom ( [4.23|) and (|5.12|) it is clear that the critical coupling Qc is the same in both the finite 
and infinite lattice cases. However, the crtitical coupling Qnc^ in the large limit, is finite 
in the infinite lattice case and zero in the case of the finite lattice. It, therefore, follows 
that the transition to the vortex phase occurs at finite temperature for infinite lattice and 
at zero temperature for the finite lattice. 



6. Conclusion 

The noncommutative X — Y model we have considered is a generalization of the 
commutative X — Y model for nonzero noncommutativity parameter 9. By mapping the 
field theory to its equivalent matrix model we have been able to construct vortex solutions 
in it. Kosterlitz Thouless transition from a phase of no vortex to a phase of vortices 
occurs when the field theory coupling or, equivalently, the matrix model coupling, reaches 
a critical value. We have shown that the critical matrix model coupling is the same when 
the system is put in either a finite or an infinite lattice. The critical values of the coupling 
in the continuum field theory is finite in the case of an infinite lattice and zero in the case 
of a finite lattice. 
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